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ABSTRACT 

We investigate the thermal and ionization history of Lya clouds photoionized by a time- 
dependent UV background, including nonequilibrium effects. The results show that it 
is possible to obtain temperatures as low as T ~ 15000 K (or, equivalently, Doppler 
parameters b ~ 15 km s _1 ) at z = 3 for cloud total densities n ~ 10~ 4 cm -3 , if (i) 
the reionization epoch occurred at z% ~ 10, and (ii) the UV background has a factor 
70—100 decrease at the Hell edge. A trend towards smaller b with increasing redshift is 
present in the redshift interval 2 = 1 — 5. Higher densities lead to higher values of b and 
smaller hydrogen correction factors, nnii/nm< The correction factors for helium are 
also given. For a hydrogen column density Nhi = 3 x 10 14 cm~ 2 , cloud sizes are larger 
than 100 kpc, consistent with recent observations of quasar pairs. Pressure confined 
models, instead, yield implausibly low cloud densities at low redshift, and too small 
sizes at intermediate redshift. The implications of the model are confronted with the 
available observational data. 



1 INTRODUCTION 

The study of the Lya absorbers seen in quasar spectra short- 
ward of the quasar Lya emission has a considerable cosmo- 
logical relevance as a probe of the physical properties of the 
neutral gas phase of the baryonic component of the Universe. 

The optically thin Lya lines are usually associated to 
highly ionized clouds (un/uni ~ 10 4 ) in photoionization 
equilibrium with the general UV ionizing background (UVB) 
produced by quasars. The clouds may reside in the inter- 
galactic space between galaxies or be associated with the big 
halos surrounding high redshift galaxies (Bahcall & Spitzer 
1969; Sargent et al. 1980; Lanzetta et al. 1994). 

The knowledge of the physical and cosmological prop- 
erties of the Lya absorbers relies on a few parameters (red- 
shift, column density, Doppler width) which can be derived 
from the line-profile fitting of (often blended) absorption 
features. 

The column density distribution is a double power-law 
with a break at log Nhi ~ 14 below which the slope is /3f ~ 
1.4. Above the break the slope is steeper, /3 S ~ 1.8 (Giallongo 
et al. 1996). 

High resolution data taken at resolution R ~ 30000 
seem to indicate typical Doppler values of the order of 
b ~ 23 — 25 km s _1 corresponding to cloud temperatures 
T ~ 32000 - 38000 K (Cristiani et al. 1995, Hu et al. 1995, 
Giallongo et al. 1996). 



It is important to note that narrow absorption lines with 
b ^ 20 km s -1 have been observed in high resolution Lya 
samples. Chaffee et al. (1983) first noted that a line in their 
spectrum at resolution ~ 12 km s _1 was surprisingly narrow 
with b ~ 15 km s _1 . More recent echelle data with variable 
signal-to-noise ratios reveal a non negligible fraction of nar- 
row lines with Doppler parameter in the range 10 < b < 20 
km s _1 . The few lines (~ 1%) with b ^ 10 km s _1 present 
in the spectra have been often identified as metal-lines asso- 
ciated with heavy element systems (e.g. Rauch et al. 1993, 
Hu et al. 1995). 

It is difficult to derive the intrinsic fraction of narrow 
lines with b ~ 15 — 20 km s _1 since it depends on the 
S/N ratio, intrinsic blending of the lines, and resolution of 
the various spectra and on the fitting procedure adopted. 
These fractions range from 10% to 20% of the lines with 
log Nhi Z 12.8 (Carswell et al. 1991; Rauch et al. 1993; 
Cristiani et al. 1995; Hu et al. 1995; Tytler 1995; Giallongo 
et al. 1996). 

In any case, reconciling temperatures as low as T ~ 
15000 K with sizes along the l.o.s. greater than 30-50 kpc 
for the Lya clouds proves to be a difficult challenge. If the 
clouds are assumed in photoionization and thermal equilib- 
rium with a power-law UV background produced by quasars 
at z = 2 — 4, then temperatures T < 30000 K could 
only be obtained increasing the hydrogen volume density 
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to uh > 10 3 cm 3 . As a consequence, cloud sizes along 
the line of sight drop down to Su < 100 pc, as outlined by 
Chaffee et al. (1983), Pettini et al. (1990), and Donahue & 
Shull (1991). 

Besides, a growing number of Lya forest observations 
in quasar pairs suggest large transverse sizes of the clouds. 
Foltz et al. (1984) derived a lower limit in the range 5 — 
25 ZigQ 1 kpc for the gravitational lensed pair Q2345+0007, 
while Smette et al. (1992) gave lower limits in the range 
12 - 50 h^o kpc for the QSO pair UM673. Recent measures 
by Bechtold et al. (1994) give larger values with a 90% lower 
limit 80/i^q 1 kpc at z = 1.8 for the diameters of spherical 
clouds. Smette etal (1995) find a 2a lower limit of lOOh^ 
kpc for the diameter of spherical clouds at z ~ 2. 

Given the implied low densities, radiative recombina- 
tion can not be the main cooling agent for clouds with 
T < 20000 K. Duncan, Vishniac & Ostriker (1991) proposed 
adiabatic cooling due to cloud expansion against the exter- 
nal pressure of a diffuse but hot (T > 10 J K) intergalactic 
medium. They obtained cloud temperatures T ~ 15000 K 
with sizes ~ 20 kpc. However, their calculations are based on 
the assumption of ionization equilibrium, and neglect heat- 
ing from helium photoionization, which increases the cloud 
temperature by a factor 1.8. 

Giallongo & Petitjean (1994) (GP) stressed the impor- 
tance of the inverse Compton cooling on the CMBR, and of 
a steepening of the ionizing UV background at wavelengths 
shorter than the Hell edge. The combination of the two ef- 
fects was shown to produce temperatures T ~ 22000 K with 
sizes up to 100 - 200 kpc. 

However, GP imposed both thermal and ionization 
equilibrium for clouds with densities uh ~ 10~ 4 cm -3 at 
z ~ 3, an assumption which is only marginally consistent 
with the relevant timescales of the problem. 

Models which include the effect of gravity have also been 
proposed (Petitjean et al. 1994; Meiksin 1994) in the context 
of the so called minihalo model (Rees 1986). Typical Doppler 
values > 20 km s _1 are found together with cloud sizes < 60 
kpc. 

In the last period, several authors have studied the 
properties of the Lya clouds in the general cosmological cold 
dark matter scenario (Hernquist et al. 1995; Zhang et al. 
1995; Miralda-Escude et al. 1995). 

Although these models follow the formation and evolu- 
tion of the Lya clouds in a dynamical scenario, the treat- 
ment of photoionization is only approximate either because 
of the assumption of ionization equilibrium or because of 
the uniform ionizing background adopted with fixed initial 
conditions. 

It is important to investigate time dependent photoion- 
ization models for the Lya clouds for the following reasons: 
(i) the UVB produced by quasars has a strong dependence 
on redshift; (ii) the inverse Compton cooling time, t c , is 
longer than the Hubble time, tn, at the typical redshifts 
at which Lya clouds are observed: t c /tH ~ 40(1 + z)~ 2 ' 5 
(GP). 

The aim of this study is to investigate in detail the ther- 
mal and ionization history of the clouds through a nonequi- 
librium photoionization model. This simple model allows us 
to explore the ionization state of the cloud as a function of a 
time-dependent UV background with different reionization 



epochs and with different shapes in frequency near the Hell 
edge. 

Our results show that it is possible to obtain temper- 
atures as low as T ~ 15000 K, for cloud total densities 
n ~ 10~ 4 cm~ 3 if the reionization epoch occurred at z% 10 
and if the UV background has a factor 70 — 100 decrease at 
the Hell edge. 



2 PHOTOIONIZATION MODEL 

The equations governing the thermal and ionization evolu- 
tion of a Lya cloud immersed in an isotropic ionizing back- 
ground of intensity I v are (see Ferrara & Field 1994 for de- 
tails) : 



^n X i+i = n X i-y p (X l )+'y c (X l ,T)n xt n e -n X i+in e a(X l ,T), 
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The gas is assumed of primordial composition with a 
helium abundance equal to riHe = O.Iuh and total number 
density n; even if metals are present, they do not affect the 
temperature very much provided the ionization level is high, 
as the one found here. 

We consider the evolution of the following species: 
H°,H + ,He°,He + ,He ++ , and electrons. The symbol X = 
H, He denotes the element considered, and i its state of 
ionization; obviously, i = for hydrogen, and i — 0, 1 for 
helium; n e is the electron density, 7 is the specific heat ratio 
assumed to be constant and equal to 5/3. 

The photoionization rate y p is given by 

f°° T 

J"L X W 

J v is the first moment of the field, vl x indicates the ion- 
ization limit for each species, a is the photoionization cross- 
section, (j> is the secondary ionization rate, 7 C is the colli- 
sional ionization coefficient, and a is the total recombination 
coefficient. We have adopted the "on the spot" approxima- 
tion in which the diffuse field photons are supposed to be 
absorbed close to the point where they have been generated. 

The energy equation (2.2) governs the evolution of the 
gas temperature T. We have used the following notation: 
hx = nx /n is the relative abundance of the element X; 
xx,i+i = tl^ 1 jnx is the fractional density of the ionization 
state i + 1 of the element X; fix = {imax + if [l + Y^^" (i + 
l)sx,i+i] _1 , with imax = 0, 1 for H and He, respectively. 

The function C represents the net cooling rate per unit 
volume and is given by the difference between photoioniza- 
tion heating and radiative cooling: 



Jv_ 
I IV 



The following processes have been included in the cal- 
culation of the cooling function: (i) free-free from all ions; 
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(ii) H and He recombinations; (iii) electron impact ioniza- 
tion of H and He; (iv) electron impact excitation of H and 
He (n=2,3,4 triplets); (v) He dielectronic recombination. At 
high temperatures or early epochs inverse Compton cool- 
ing on the CMWB photons is important and it has been 
included as in Ikeuchi & Ostriker (1986). 

We mostly concentrate on constant cloud total density 
cases, (i.e., dn/dt = in eq. [2.2]); in addition, we also ex- 
plore cases in which Lya clouds are taken to be in pressure 
equilibrium with (and hence confined by) the surrounding 
IGM. To this aim, it is useful to write eq. (2.2) as follows: 



(7- 



_ h d_ 

1) dt 
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n 
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h x dpi 



HX,; 



dt ' 
(2.3) 

where pi(z) is the pressure of the IGM, which is a pre- 
scribed function of the redshift z. Ikeuchi & Ostriker (1986) 
(and more recently Shapiro 1995) have studied in detail the 
thermal history of the IGM from the era of the reheating 
to the present; they conclude that the most plausible sce- 
nario requires that both photoionization and bulk mechani- 
cal heating (i.e. shocks) contribute to the heating. Ikeuchi & 
Ostriker show that, if the reionization occurred at z, = 10 - 
for example - the initial evolution of the IGM is practically 
isothermal, pi oc (1 + z) 3 , and the IGM enters the adiabatic 
expansion phase, pi oc (1 + z) 5 , at z ^ 4. As for the local 
value of the pressure, we take p/(0) = 10" 2 cm" 3 K; we will 
use the cosmology Q — 1 and hso = 1- 

The problem is completely specified once the evolution 
of the UVB field is assigned together with the initial con- 
ditions (the latter are discussed below). We have assumed 
an evolving UVB following Madau (1991) and Meiksin & 
Madau (1993) model. The intensity of the UVB at the Ly- 
man limit, Jll, evolves as (1 + z) 3 for z < 2; it is approx- 
imately constant, i.e. oc (1 + z) 1 / 4 , in the range 2 < z < 3, 
and it declines as oc exp[— 0.69(z — 3)] for z > 3. 

We assume a normalization value Jll — 7 x 10 -22 erg 
s" 1 cm -2 sr" 1 Hz -1 at z = 2, consistent with the average 
value Jll — 5 x 1CP 22 found from the proximity effect of a 
Lya high resolution sample evaluated in the redshift range 
z = 1.7 — 4.1 (Giallongo et al. 1996). A power-law spectrum 
with index a = —1.5 for the UVB has been adopted. 

In some cases a factor 100 break has been introduced at 
4 Ryd to simulate the Hell absorption trough, consistently 
with the detections of the Hell Lya Gunn-Peterson effect 
which have provided an indirect evidence of a steepening at 
the Hell edge with a UVB ratio between the Lyman limit 
and the Hell edge S L ~ 100 (Jakobsen et al. 1994). 



Figure 1. (a) Evolution of the Dopplcr perameter b (lower panel) 
and of the hydrogen correction factor Xh.i (upper panel) for con- 
stant cloud density n = 10" 4 cm" 3 as a function of redshift for 
different initial conditions (see text). The initial conditions are 
Zi = 7, 10 and Tj = 10 4 , 10 s K. Points refer to b obtained assum- 
ing ionization and thermal equilibrium for comparison assuming 
the same UVB. Solid lines indicate cases with B = 1 (no break); 
dotted lines indicate cases with B = 100. The corresponding tem- 
perature scale in the lower panel is shown on the right axis, (b) 
As Fig. la, with initial conditions z; = 7 and T; = 10 4 K, B = 1 
for different values of the density: n = 10" 3 cm" 3 (dotted line), 



10" 



cm 3 (solid), n = 10 4,5 cm 3 (long-dashed). 



3 MODEL RESULTS 

Using the photoionization model described in the previous 
section we have first explored the thermal and ionization his- 
tory of clouds having constant total density. This isochoric 
scenario, although unrealistic if followed for the overall Hub- 
ble time, allows a simple insight on the evolution of the ion- 
ization and thermal state of low density clouds in response 
to changes of the UVB. 

Among the several cases investigated, we present the 
most relevant ones featuring (i) different initial conditions 



(reionization epoch, z» = 7, 10, and temperature, 
log Ti = 4, 5) and (ii) different values of the Hell break at 4 
Ryd, B = 1, 100, assuming a total density n = 10~ 4 cm" 3 . 
The gas is supposed to be mostly neutral at z = z%; different 
assumptions produce slight changes only for redshifts out- 
side the currently observable range. At z = 10 the Compton 
cooling time is so short that there are no relevant differences 
between the two different T». 

Fig. la shows the thermal and ionization evolution 
of a Lya cloud as described by the Doppler parameter, 
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b = (2fcT/mff) 1//2 and the hydrogen correction factor, 
Xh,i = riHii/riHi = x H ,i/(l - x H ,i)- Tne Doppler param- 
eter evolves strongly as a function of redshift and, indepen- 
dently of Tt and z;, reaches a maximum value b ~ 28 km s _1 
at about z — 1.2 for B = 1, i.e., no Hell break; the memory 
of the initial conditions is completely washed away already 
at z ~ 2, for both cases B = 1 and £> = 100. The tem- 
perature peak is due to the lower efficiency of the Compton 
cooling which decreases as (1 + z) 4 ; for z 1.2 this effect 
is overwhelmed by the decrease of the UVB, and hence of 
the photoionization heating. The peak position also lags be- 
hind an analogous maximum of xh,i (see below) since the 
decrease in the photoionization heating is initially compen- 
sated by a larger neutral fraction. 

The cases B = 100 look qualitatively similar to the 
B = 1 ones, but they show a general tendency towards lower 
temperatures, with a minimum value b ~ 13 km s _1 attained 
at z = 6 — 7, with z, = 10. It is worth mentioning, though, 
that even at z = 3 typical b values are 16 — 18 km s _1 , to 
be compared with b ~ 23 km s _1 for the no break case. 

The hydrogen correction factor Xha traces to a good 
extent the evolution of the UVB, peaking at z = 2; the range 
of variation of Xh,i is 3000 — 10 5 . The evolution is rather 
insensitive to z, and Tj, the different curves converging as 
early as z ~ 5. A B — 100 break reduces Xh,i by ~ 35% 
around the maximum (z = 2). 

We have also explored the cases n — 10~ 3 cm~ 3 and 
n = 10~ 4,5 cm~ 3 with initial conditions z, = 7 and B = 1; 
the corresponding evolutionary curves are shown in Fig. lb. 
The behavior Xh,i is qualitatively the same as above; in par- 
ticular the value of the maximum increases almost linearly 
with n. The Doppler parameters are lower that for the stan- 
dard n = 10~ 4 cm~ 3 case both for n — 10 -3 cm~ 3 and for 
n — 10~ 4,5 cm~ 3 (this is due to a transition from Compton 
cooling dominated losses at low n to recombination dom- 
inated cooling at higher n). In addition, the maximum of 
the b curve shifts towards lower redshifts as the density is 
decreased. 

In general, we can define correction factors for helium 
in analogy to the hydrogen one as XHe,i+i = riHe,i+i/nHi = 

h, HeXHe,i+i/hH(l — xh,i) (notice the definition relative to 
neutral hydrogen). Fig. 2 shows the behavior of the helium 
correction factors as a function of redshift. The most relevant 
feature is the change in the evolution of the Hell fraction 
from the cases B = 1 to the cases B — 100. For B = 1 
and z < 5, XHe,i has an almost constant low value ~ 20, 

i. e., most of the He is in the doubly ionized state. As B is 
increased to 100, Xffe.i develops a maximum around z = 3 
of about 4000, a value similar to that of Xh&,i in the same 
redshift interval. Thus the ratio of the two correction factors 
provides a very sensitive measure (and test) of the 4 Ryd 
break. 

The previous results allow to compute the Lya cloud 
sizes along the l.o.s. Sy = Nhi{1 + XH,i)/n at any given 
redshift, as illustrated by Fig. 3. Typical sizes of the clouds 
are several kpc, ranging from a minimum of about 5 kpc for 
N H i = 10 13 cm" 2 to 100-120 kpc for A H / =3xl0 14 cm" 2 in 
the redshift interval 2-3. Given the higher neutral fraction 
found in models with B = 100, sizes are a factor ~ 1.5 
smaller. Apart from B, the sizes are quite insensitive to the 
initial conditions. Locally (z = 0), clouds with Nhi < 3 x 
10 14 cm -2 are predicted to have sizes not exceeding 5 kpc. 



Figure 2. Evolution of the helium correction factors, XHe.i an d 
XHe.,2 for constant cloud density n = 10~ 4 cm -3 as a function of 
redshift for the same cases as in Fig. la. Solid lines indicate cases 
with B = 1 (no break); dotted lines indicate cases with B = 100. 

Finally, models in which Lya clouds are assumed to 
be in pressure equilibrium with the IGM, according to eq. 
(2.3), are presented in Fig. 4. The temperature of the clouds 
is now a monotonically increasing function of time, as the 
cooling reduction caused by the density decrease dominates 
the opposite effect due to the decrease of the photoioniza- 
tion energy input. At z = 3, b = 23 — 25 km s _1 , depend- 
ing on B. We point out that these results depend on the 
particular assumption made for the thermal history of the 
IGM, following Ostriker & Ikeuchi (1986); deviations from 
these values may produce quantitative changes in the deter- 
mination of b. The same behavior is shown by Xh,i, which 
reaches very high values (~ 3 x 10 6 ) at z = 0. The total den- 
sity n of the cloud drops very rapidly from the initial value 
rii ~ 10~ 2 cm~ 3 to an extremely low value ~ 10~ 7 cm~ 3 
at z = (essentially the IGM density). As a consequence, 
sizes are implausibly large, or, stated in a different way, the 
clouds would dissolve rapidly below z ~ 2. At intermedi- 
ate redshifts, pressure confined models yield essentially the 
same b values as the constant density, B = 1 cases (but 
higher than B = 100 ones) ; 

however, correction factors (and thus sizes) are typically 
one order of magnitude smaller. Thus, low-z Lya clouds ob- 
servations can be used to put strong constraints on pressure 
confinement models. 

4 DISCUSSION AND IMPLICATIONS 

We first compare our time-dependent results with the ones 
obtained by GP who assumed thermal and ionization equi- 



PROPERTIES OF THE LYMAN ALPHA CLOUDS 5 



Figure 3. Cloud sizes Sy for constant cloud density n = 
10~ 4 cm~ 3 as a function of redshift for initial conditions = 10, 
Ti = 10 4 K. Solid lines correspond to B = 1 (no break); dotted 
lines correspond to B = 100. Three different hydrogen column 
densities Nhi = 10 13 ,10 14 ,3 X 10 14 cm~ 2 are shown for each 
case. 

librium. GP used a value of Jll = 1CF 21 erg s _1 cm~ 2 
Hz -1 , a value corresponding to z = 2.25 for the presently 
adopted J v {z). At that redshift, the cloud temperature (see 
Fig. la) is T ~ 3 x 10 4 K for B = 1, and T ~ 2 x 10 4 K 
for £> = 100, to be compared with the analogous values ob- 
tained by GP (see their Fig. 1), T ~ 4.7 x 10 4 K for B = 1, 
and T ~ 3 x 10 4 K for B = 100. We have plotted in Fig. 
la the values of b obtained assuming ionization and thermal 
equilibrium for comparison assuming the same UVB. The 
reasons for the discrepancy are: (i) Compton losses on the 
CMWB strongly cool the clouds at high redshift; (ii) the 
UVB decreases exponentially for z > 3; and (iii) nonequilib- 
rium effects. The sizes found by GP are about a factor 1.5 
larger for the same Nhi and n, thus implying a larger hydro- 
gen correction factor. This is due to the enhanced collisional 
ionization produced by their higher temperature. 

The bottom line of our study is that it is feasible to 
reconcile two apparently conflicting observed properties of 
the Lya clouds, low temperatures and large sizes, basi- 
cally with one assumption, i.e. the existence of a break 
in the UVB at the Hell edge. It is easy to show that 
Lya clouds become optically thick to radiation above the 
Hell Lyman limit for a neutral hydrogen column density 
Nhi = 6.3 x 10 17 / X He,i cm . For the calculated values of 
XHe,i ~ 20 (B = 1), Lya clouds with Nhi > 3 x 10 16 cm~ 2 
will substantially absorbe the flux above 4 Ryd, thus enhanc- 
ing the value of XHe,i even further (see Fig. 2 for B = 100). 
This well-known (Miralda-Escude & Ostriker 1990, Shapiro 
1995) nonlinear effect is thus responsible not only for the 



Figure 4. As Fig. la, but for the constant pressure case, and 
initial conditions Zj = 7, 10, Tj = 10 4 K. The dashed line in 
the upper panel is the total density, n, of the cloud. Solid lines 
correspond to B = 1 (no break); dotted lines correspond to B = 
100. 



break, but also, given the reduced Hell ionizing flux, for a 
detectable Hell Gunn- Peterson effect (Jakobsen et al. 1994). 
In addition, from our results we do expect the presence of a 
pronounced He II Lya forest. 

The issue about Doppler parameters and cloud sizes 
deserves some additional discussion. First it appears that b 
values lower than about 15 km s _1 cannot be reproduced by 
any of our nonequilibrium models. This agrees with the pres- 
ence of a lower envelope of the 6-distribution with a cut-off 
at b c ~ 16 km s _1 found by different investigators (Car- 
swell et al. 1991, Rauch et al. 1992, 1993, Hu et al. 1995, 
Giallongo et al. 1996), disputing a previous claim by Pettini 
et al. (1990) of a lower value of b c , probably affected by line 
selection and fitting procedure artifacts. Quite interestingly, 
we find also an upper limit b ~ 28 km s _1 , whereas the 
above studies have inferred b- values as high as 60 km s _1 . 
The most natural suggestion is that high-6 values can be un- 
resolved blends of b ~ b c lines (Rauch et al. 1992, Hu et al. 
1995). However, in this case the two-point correlation func- 
tion should display an increased amplitude on small scales, 
typically Av < 100 km s _1 . Some indirect support to this 
hypothesis is lent by Cristiani et al. 1995 and Meiksin & 
Bouchet (1995), who both found significant clustering on 
scales Av > 100 km s _1 . If blending is viable, depatures 
from Maxwellian velocity distribution are expected (Kulka- 
rni & Fall 1995). A different explanation, i.e. bulk motions, 
has been shown (Press & Rybicki 1993) to encounter ap- 
parently insurmountable difficulties. Alternatively, if a hot 
IGM does exist, temperatures of the order of a few x 10 5 K 
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could be generated in the interfaces between the IGM and 
the clouds by thermal conduction (Ferrara & Shchekinov 
1996). 

The cloud sizes S\\ derived here for n = 1CF 4 cm~ 3 , al- 
beit larger than in previous models, are still smaller than 
the transverse ones, S±, inferred either from lensed or 
paired quasars. Thus, unless cloud densities are smaller than 
~ 1CT 4 cm -3 , we are forced to conclude that clouds must 
be flattened, a conclusion shared by other investigators on 
different grounds (Rauch & Haehnelt 1995, Hu et al. 1995). 
Dinshaw et al. (1994,1995) have provided very interesting 
evidence for the evolution of S±: from their quasar pair ex- 
periments they obtain S± > 320/igQ 1 kpc in z = 0.5 — 0.9 
and S± > 80/ig, 1 kpc at z — 1.8. From Fig. 3, independently 
of N H i, the ratio S\\(z = 1.8)/S||(z = 0.5) = 5.6, while the 
analogous ratio for S± is ~ 0.25. Although very speculative, 
this finding seems to suggest that the clouds evolve from an 
approximately round shape to a highly flattened one, remi- 
niscent of a gravitational collapse, between z ~ 2 and z ~ 1. 
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